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1.  INTRODUCTION 


A  collaborative  program  to  test  the  hypothesis  that  electromagnetic  fluctuations  at  the  sea 
floor  are  forced  by  turbulent  vorticity  fluctuations  in  the  bottom  boundary  layer  began  in  Fall, 
1993.  After  testing  of  the  instrumentation  in  the  early  part  of  the  year,  vorticity  data  were 
recorded,  together  with  electromagnetic  fluctuations,  in  the  strong  tidal  flows  in  the  Strait  of  Juan 
de  Fuca  through  four  seasons.  Summer  1994  through  Spring  1995.  The  vorticity  measurements 
provide  data  for  testing  existing  theoretical  predictions  of  the  spectrum  of  the  vorticity 
fluctuations  at  wave  numbers  within  the  inertial  subrange,  and  they  can  provide  a  means  to  charac¬ 
terize  the  vorticity  spectrum  at  wave  numbers  below  the  inertial  subrange. 

The  Woods  Hole  Oceanographic  Institution  (WHOI)  component  of  the  two-year  program 
was  to  provide  measurements  of  vorticity  at  two  scales,  45  cm  and  150  cm,  together  with  measure¬ 
ments  of  horizontal  velocities  between  1  and  3.4  m  above  bottom.  The  Scripps  Oceanographic 
Institution  (SIO)  component  of  the  program  was  to  simultaneously  measure  electromagnetic  fluc¬ 
tuations  near  the  WHOI  vorticity  meters.  SIO  also  observed  velocity  from  three  nearby  moorings 
between  3  m  and  30  m  above  bottom.  Conductivity-temperature-depth  (CTD)  casts  were  taken 
by  SIO  during  deployment  and  several  months  later.  A  side-scan  sonar  survey  of  the  region  within 
a  tidal  excursion  of  the  deployment  site  was  conducted  by  WHOI  in  July,  1995.  The  program  was 
funded  by  the  U.  S.  Office  of  Naval  Research  and  was  aimed  at  providing  technology  which  can 
enhance  the  performance  of  existing  nonacoustic  surveillance  instrumentation  with  the  ability  to 
incorporate  environmental  knowledge  at  the  observation  location. 

2.  INSTRUMENTATION 


Chip  Cox  and  Dave  Jacobs,  of  SIO,  designed  and  built  an  instrument  to  observe 
electromagnetic  fluctuations.  The  instrument,  named  Charlotte’s  Webb,  has  six  radial  arms  with 
pressure  and  electric-field  sensors  at  0, 1.2, 2.4, 4.8  and  9.8  meters  fi-om  the  center.  There  are  also 
sensors  on  the  midpoint  of  the  chord  between  spokes  at  4.8  m  radius.  Two  pair  of  electrodes  were 
mounted  vertically,  1  m  above  the  center  point.  SIO  also  deployed  three  moored  current  meters 
and  conducted  CTD  casts  during  deployment  and  recovery  of  a  small  "modem  pod". 

WHOI  designed  and  built  two  tripods  which  recorded  temperature,  pressure,  velocity  and 
vorticity.  The  tripods  were  also  equipped  to  measure  orientation,  pitch  and  roll  of  each  tripod.  The 
tripods  are  illustrated  in  Figures  1  and  2  and  described  in  Tables  1  and  2.  The  vorticity  meters  were 
designed  using  electronics  based  on  the  Benthic  Acoustic  Stress  Sensor  (BASS)  (Williams  et  al., 
1987). 


As  illustrated  in  Figure  1,  VORT#2  housed  two  small-scale  vorticity  meters  (45  cm)  to 
provide  horizontal  components  of  velocity  at  1.0, 1.4,  3.0  and  3.4  m  above  bottom,  as  well  as 
three-dimensional  vorticity  at  1.2  m  and  3.2  m  above  bottom.  VORT#3,  as  illustrated  in  Figure 
2,  housed  one  large-scale  vorticity  meter  (150  cm),  which  was  built  within  the  structure  of  the 
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tripod,  providing  horizontal  components  of  velocity  at  1.1  m  and  2.4  m  above  bottom  and  vorticity 
at  1.25  m.  For  each  vorticity  meter,  the  velocity  was  observed  along  twelve  acoustic  paths,  defined 
as  ul,  u2,  u3  ...,  ul2,  as  shown  in  Figure  3a.  These  paths  can  be  used  to  compute  horizontal 
velocity  (Figure  3b)  and  vorticity  (Figure  3c).  The  use  of  each  path  and  the  definition  of  a 
real-world  coordinate  system  are  described  in  Section  5. 

The  differential  pressure  sensors  were  Paroscientific^  Model  No.  8DP700-1.  The 
manufacturer’s  calibration  coefficients  and  algorithm  are  included  in  Appendix  A  as  seepress.m. 
On  VORT#2,  thermistors  were  potted  into  the  vorticity  meter  stalks,  about  2.5  cm  away  firom  the 
transducer.  VORT#3  thermistors  were  strapped  to  a  leg.  The  thermistors  (YSI^)  were  calibrated 
at  WHOI,  as  described  in  Section  5. 

Each  tripod  had  been  equipped  with  an  Ocean  Instrument  Systems^  acoustic  release  mech¬ 
anism  and  a  Benthos^  sphere  ("pop-up")  housing  a  hard  drive  which  recorded  mirrored  data  from 
each  data  logger  using  inductive  modems.  Each  sphere  was  also  equipped  with  an  ARGOS^ 
transmitter  to  assist  in  locating  the  pop-up  upon  release. 


1. Paroscientific,  Inc.  Redmond,  WA  98052 

2.  Yellow  Spring  Instruments. 

3, Ocean  Instrument  Systems,  North  Falmouth,  MA  02556 

4. Benthos,  Inc.,  North  Falmouth,  MA  02556 

5,  ARGOS  receiver/transmitters  are  manufactured  by  NACLS,  INC.,  Landover,  MD  20785 
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HEIGHT 
OFF  BOTTOM 


Figure  1.  VORT#2  small  scale  vorticity  meter,  showing  the  heights  of  the  horizon 
tal  velocity  planes.  The  length  of  each  acoustic  path  is  45. 


Figure  2.  VORT#3  large  scale  vorticity  meter,  showing  the  heights  of  the  horizontal 
velocity  planes.  The  length  of  each  acoustic  path  is  150  cm. 
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Table  1.  VORT#2  Instrumentation 

Measurement 

Height  Above 
Bottom  (cm) 

Description 

Thermistor#1  104 


Velocity 

105 

Bottom  Sensor  -  Bottom  Plane 
(Paths  1,5,9) 

Vorticity 

(Center  of  Each  Plane) 

Velocity 

■ 

Bottom  Sensor  -  Top  Plane 
(Paths  3,7,11) 

Thermistor#2  149 

Pressure  227  Paroscientific  SN57216 

Thermistor#3  304 

Velocity  308  Top  Sensor  -  Bottom  Plane 

_ (Paths  13,17,21) _ 

Vorticity  337  (Center  of  Each  Plane) 

Velocity  350  Top  Sensor -Top  Plane 

_ (Paths  15,19,23) _ 

Thermistor#4  354 


Table  2.  VORT#3  Instrumentation 

Measurement 

Height  Above 
Bottom  (cm) 

Description 

Pressure 

33 

Paroscientific  Inst  SN57222 

Thermistor#1  36 

Thermistor#2  59 


Velocity  112  Bottom-Plane 

_ (Paths  1,5,9) _ 

Thermistor#3  112 

Thermistor#4  139 

Thermistor#5  167 

Vorticity  175  (Center  of  Each  Plane) 

Thermistor#6  210 

Velocity  240  Top-Plane  (Paths  3,7,11) 

Thermistor#7  251 

Thermistor#8  297 
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Figure  3a.  Definition  of  each  acoustic  path  for  VORT#3.  The  explanation  of  u  and  v  can  be  found  in 
Section  5. 
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Figures  3b  &  3c 
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Figure  4.  Bathymetric  map  shows  depth  contours  (in  meters)  with  location  of 
pre-deployment  test  sites.  The  tidal  excursion  was  estimated  to  be  approximately  5 
km,  based  on  tidal  velocities  of  0.75  m/s.  VORT#2  and  VORT#3  were  positioned 
approximately  60  meters  apart,  along  the  10  meter  isobath. 


3.  PRE-DEPLOYMENT  TESTS 

WHOI  conducted  three  deployments  to  test  the  newly  developed  tripods  and 
instrumentation.  The  tripods  were  placed  approximately  200  feet  apart  about  a  quarter-mile  east 
of  Nobska  Point  (Woods  Hole,  MA)  in  Vineyard  Sound  along  the  10  meter  isobath.  (See  Figure 

4. )  The  tests  proceeded  as  follows: 

May  7-10  1994  -  to  observe  instruments  for  strumming  in 
high  currents  (no  electronics) 

May  28-30, 1994  -  with  electronics  continuously  sampling  at  4  Hz 
for  less  than  one  day 

June  17-22, 1994  -  with  electronics  sampling  at  7.14  Hz  for  72  seconds, 
every  other  minute  for  34  minutes  every  third  hour 
for  five  days 

The  May  7th  deployment  indicated  that  no  significant  strumming  occurred,  as  observed  by 
three  divers  during  peak  tide,  which  is  presumed  to  be  approximately  1  m/s  (based  on  subsequent 
measurements). 

During  the  May  28th  deployment,  continuous  data  were  collected  until  the  disk  was  full, 
which  occurred  after  approximately  14  hours.  It  appears  that  the  VORT#2  tripod  tipped  over 
during  strong  flow.  Therefore,  there  were  less  than  2  hours  of  data  for  the  small  scale  vorticity 
meter.  Upon  analysis  of  the  pressure  data  from  VORT#3,  we  learned  that  a  clock  interrupt  was 
not  being  serviced,  effecting  a  clock  cycle  of  280  ms,  not  250  ms,  as  planned.  This  affects  the 
pressure  measurement  and  the  frequency  spectrum  (i.e.,  3.57  Hz  sampling,  not  4Hz).  The 
VORT#3  electronics  had  been  set  to  ±0.75  m/s  for  maximum  count  and,  therefore,  velocity  buffers 
overflowed  at  peak  current.  These  problems  have  subsequently  been  corrected.  The  internal 
clocks  started  when  power  was  supplied,  causing  a  0.5165  day  lead  in  the  VORT#2  clock. 


The  thermistors  were  calibrated  at  the  WHOI  CTD  Calibration  Lab  on  June  9-10. 
Coefficients  for  conversion  to  temperature  are  listed  in  Appendix  A,  with  the  MATLAB®^  script 
used  for  the  conversion.  The  sampling  scheme  was  similar  to  that  actually  used  in  the  Strait  of 
Juan  De  Fuca.  VORT#2  and  VORT#3  dataloggers  were  turned  on  June  17th  at  1 1 :30  and  1 1 :45 
EDT,  respectively.  The  tripods  were  deployed  at  12:16  and  12:18  EDT,  respectively.  No  data 
from  the  VORT#2  tripod  were  recorded,  due  to  an  electronics  failure.  The  acoustic  paths  on 
VORT#3  were  atypically  noisy.  This  noise  is  reasonably  explained  as  being  attributable  to 
seaweed.  The  tripod  was  covered  with  it  upon  recovery  and  there  appears  to  be  more  noise  in  the 
top-most  plane,  where  the  acoustic  paths  are  below  tubing  which  trapped  seaweed. 


6.MATLAB  is  a  product  of  The  Mathworks,  Inc.,  Natick,  MA  01760 
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Figure  5.  Bathymetric  map  shows  depth  contours  (in  meters)  with  location  of  deployment  sites.  The 
tidal  excursion  was  estimated  to  be  approximately  5  km,  based  on  tidal  velocities  of  0.75  m/s. 
VORT#2  and  VORT#3  were  positioned  approximately  600  meters  apart,  along  the  125  meter  isobath. 
The  positive  along-strait  velocity  (u)  is  defined  as  flow  along  the  isobath  (positive  is  eastward).  The 
cross-strait  velocity  (v)  is  defined  as  flow  shore-normal  flow  (positive  is  northward). 


Table  3.  Instrument  Summary 

Map  Site  / 
Organization 

Instrument 

Name 

Property 

Observed 

Latitude 

(North) 

Longitude 

(West) 

Deployment 

(GMT) 

1/SIO 

Charlotte’s 

Webb 

Electro¬ 

magnetic 

Fluctuations 

48°  23.59’ 

^24^  32.19’ 

-7/11/94 

15:00 

2/SIO 

M1 

CTD/VMCM 

48^23.56’ 

124^  32.49’ 

3/SIO 

M2 

CTD/VMCM 

48^  23.79’ 

124®  32.24’ 

4/SIO 

M3 

CTD/VMCM 

48°  23.51’ 

124°  31 .84’ 

5/WHOI 

VORT#2 

Vorticity 

Temperature 

Pressure 

48°  23.55’ 

124°  31 .91’ 

7/11/94  22:48 

6/WHOI 

VORT#3 

Vorticity 

Temperature 

Pressure 

48°  23.62’ 

124°  32.43’ 

7/11/94  20:36 

10 


4.  DEPLOYMENT  &  RECOVERY 


The  tripods  were  shipped  to  the  University  of  Washington  (UW),  where  they  were 
reassembled  and  loaded  onto  the  Research  Vessel  Thomas  G.  Thompson  and  deployed  with  the 
SIO  instruments  in  the  Strait  of  Juan  De  Fuca,  as  shown  in  Figure  5. 

The  VORT#2  sensors  were  placed  in  barrels  of  stationary,  unstratified  water  to  obtain 
pre-deployment  zeros,  which  were  recorded  for  subsequent  processing.  The  VORT#3  tripod  was 
tested  from  the  ship,  while  it  was  docked  at  UW.  No  zeros  can  be  measured  for  the  VORT#3 
tripod  due  to  its  size. 

The  clocks  for  VORT#2  and  VORT#3  were  started  on  July  11, 1994,  at  18:04  and  18:06 
GMT,  respectively.  The  SIO  CTD  casts  were  conducted  and  the  SIO  current  meters  were 
deployed  after  Charlotte’s  Web  but  before  VORT#3.  VORT#2  had  been  deployed  prior  to 
VORT#3  (at  19:42  and  19:55)  but  had  to  be  recovered  when  the  recovery  floats  released  from  the 
strain  of  being  dragged  on  deployment.  After  a  stronger  bolt  secured  the  release  mechanism, 
VORT#2  was  redeployed  at  the  position  and  time  listed  in  Table  3.  The  depth  at  the  site  was  noted 
as  125  meters,  from  the  ship’s  fathometer.  The  GPS  positions  noted  in  Table  3  have  inherent 
inaccuracies  on  the  order  of  50-100  meters. 

On  August  23,  the  pop-ups  were  triggered  and  recovered  aboard  a  vessel  (Satin  Doll) 
chartered  from  Far  West  Resorts.  The  data  were  recovered  for  the  large  scale  vorticity  meter 
(VORT#3)  and  transferred  for  processing.  It  became  apparent  at  this  time  that  VORT#2  had 
failed,  although  it  had  successfully  logged  the  zeros  the  day  before  deployment. 

The  logger  was  filled  by  March  27, 1997,  and  the  instruments  were  recovered  aboard  the 
SIO  RfV  Robert  Gordon  Sproul,  on  May  28, 1997. 
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July  1 1  -  Aug.  11,1 994  only  angles  <  5  degrees 


Figure  6a  &  6b.  Comparison  of  along-strait  velocity  components  from  paths  9  (o)  and  1 1  (x)  with 
each  corresponding  reconstructed  along-strait  velocity. 
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Figure  6c  &  6d.  Comparison  of  cross-strait  velocity  components  with  along-strait  velocity 
components  in  each  velocity  plane. 


12 


5.  DATA  PROCESSING  &  ANALYSIS 

VORT#3  data  were  unpacked  and  stored  in  files  of  daily  records  with  one  variable  per  file. 
Each  continuous  time  series  consists  of  514  samples  (72  seconds  x  7.14Hz),  measured  every  other 
minute  (seventeen  times),  providing  approximately  34  minutes  of  sampling  (514  x  17  or  8738 
records):  72  seconds  on  /  48  seconds  off.  These  observations  were  repeated  every  other  hour  for 
24  hours  (104856  records)  every  third  day  for  a  month.  This  sampling  was  repeated  every  third 
month  until  the  disk  was  full,  which  occurred  in  March  1995.  The  compass  heading,  pitch  and  roll 
remained  consistent  at  332°N,  -2°  and  4°,  respectively.  The  unpacking  and  storage  routine, 
bin2mat.c,  is  included  in  Appendix  C  with  the  TTBASIC®  programs  which  controlled  the 
sampling. 

A  coordinate  system  is  defined  so  that  x  is  along-strait  (positive  in  the  flood  direction),  y  is 
cross-strait  (positive  is  poleward,  +21°Ntn,e),  and  z  is  vertical  (positive  upward).  The  correspond¬ 
ing  velocity  components  are  u,  v  and  w,  respectively.  It  can  be  shown  that 


Ubot  =  u9, 

Utop  =  -Ull, 

where  subscript  hot  and  top  refer  to  the  horizontal  planes,  as  defined  in  Table  2  and  illustrated  in 
Figure  3b,  and  uN  represents  the  acoustic  path  as  defined  in  Figure  3a. 

The  cross-strait  velocity  (v)  were  computed  as  follows: 

Vbot=  (ul-u5)/(2*cos(30°)), 

Vtop=  (u7-u3)/(2*cos(30°)). 

The  combination  of  the  other  two  horizontal  paths  in  each  plane  provides  a  consistency 

check: 


ucbot  =  -(ul+u5), 
uctop  =  (u7+u3). 

Comparisons  of  the  along-strait  velocities  (u)  with  the  reconstructed  along-strait  velocities  (uc)  are 
presented  in  Figures  6a  and  6b.  Figures  6c  and  6d  show  the  cross-strait  velocity  (v)  with  the 
along-strait  velocity.  When  the  difference  between  the  observed  velocity  and  the  reconstructed 
velocity  are  plotted  with  the  angle  of  flow  (Figure  7a  and  7b) ,  it  appears  that  there  is  a  flow 
disturbance  when  the  direction  of  flow  deviates  by  greater  than  5  degrees  from  the  u-path 
direction.  This  disturbance  is  also  apparent  in  the  estimates  of  drag  coefficient  firom  mean  velocity 
(u)  log-profiles,  as  seen  in  Figures  8a  and  8b.  For  this  reason,  we  have  restricted  some  of  the 
analyses  to  those  periods  of  undisturbed  flow.  A  more  detailed  discussion  of  flow  error  in  the 
vorticity  meter  can  be  found  in  Thwaites  (1995). 
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According  to  the  Kolmogorov  model  (Batchelor,  1953),  the  spectral  density  within  the 
inertial  subrange  can  be  used  to  estimate  the  dissipation  rate  as  follows: 

S(k)  = 

where  the  constant  P  depends  on  whether  the  spectral  estimate  is  determined  from  the  along-path 
flow  (which  is  the  case  for  paths  9  and  1 1)  or  from  the  across-path  flow  (paths  10  and  12)  and  is 
defined  as  follows: 

P  =  9/55  for  paths  9  and  11, 
and 

P  =  12/55  for  paths  10  and  12. 

As  seen  in  Figures  9-11,  filtering  of  the  turbulence  is  exhibited  when  the  flow  is  in  the 
direction  of  the  acoustic  path.  Although  this  effect  is  not  consistent  with  quantitative  calculations 
(Kaimal  et  al.,  1968),  it  appears  to  be  repeatable.  To  avoid  influence  of  the  bottom  and  contribu¬ 
tion  from  aliasing,  wave  numbers  for  the  estimate  were  restricted  to  kz  >  5  and  kz  <  kj/lO,  where 
ks  =  27tFs/u. 

and  Fs  is  the  sampling  frequency.  Therefore,  we  assumed  an  along-path  filter  proportional  to  1/k 
and  computed  an  estimate  of  a  dissipation-like  quantity  as: 

Au9,uii  =  2:(S(k)-k^'^-k)  /Nk=  S(S(k)-k^'2)/Nk 
evaluated  fromkz>5  andkz<ks/10. 


It  can  be  shown  that  the  vortex  shedding  in  the  cross-flow  paths  (ul0/ul2),  is  consistent 
with  tubing  diameters  of  the  tripod  structure  and  also  coincides  with  flow  direction,  i.e.  path  ulO 
is  contaminated  strongly  during  strong  ebb  and  ul2  is  most  affected  during  strong  flood.  For 
cross-path  flow,  ul0/ul2,  spectral  density  analysis  is  restricted  to  regions  above  kz  =  5,  and  is 
further  restricted  by  the  apparent  vortex  shedding  to  regions  below  kz  =  12.  For  the  cross-path 
estimates, 

Auio,ui2  =  i:(S(k)-k^^^)/Nk. 


As  well  as  using  only  periods  of  undisturbed  flow,  for  the  estimates  of  dissipation,  these 
data  sets  were  further  restricted  to  periods  of  stationary  flow.  The  reverse  arrangements  test,  as 
described  in  Bendat  (1986),  determined  stationarity  in  the  mean  and  standard  deviation  of  the 
mean  velocity,  u,  at  2.4  m  above  bottom.  Summaries  of  the  estimates  of  A  are  included  as  Figures 
12a  and  12b. 


14 


To  compute  vorticity,  the  four  acoustic  paths  which  form  a  closed  square  are  summed  and  divided 
by  the  length  of  each  path  (1.5  for  the  large  scale  vorticity  meter,  VORT#3): 

Va  =  (ul+u2+u3+u4)/L; 

Vb  =  (u5+u6+u7+u8)/L; 

Vc  =  (u9+ulO+ull+ul2)/L; 

The  vorticity  at  each  sample  were  rotated  into  the  x,  y,  z  coordinate  system  as  described  in  Appen¬ 
dix  B.  The  48  seconds  of  unobserved  vorticity  were  included  as  zeros  and  spectral  analysis  of  the 
vorticity  was  computed  using  4  minute  windows.  July  spectra  are  shown  in  Figure  13. 

An  estimate  of  the  variance  in  vorticity  in  a  given  observation  was  determined  by 

kz<\ 

J  5(fc)  •  2dk 

fcz>0.1 

Figure  14a  and  14b  confirm  the  nonconformity  of  ebb  and  flood  conditions  in  vorticity,  as  seen  in 
the  velocity  and  dissipation  analyses. 
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Figure  12a.  u9(bottom-o)  &  ull(top“X) 

(dashed  line  represents  a  slope  of  9/55aCd) 


Figure  12b.  ul0/ul2  (dashed  line  represents  a  slope  of  12/55aCd) 


Figure  12.  oc  =  1.5  and,  based  on  Figure  8b,  Cd(flood)  —  0.0123,  Cd(cbb)  —  0.0035U‘|u|/(kz) 


is:;,  Sv{k)*2dk 


6.  DATA  SUMMARIES 


Figures  15-24  provide  time  series  of  the  horizontal  components  of  velocity,  along  with  the 
temperature  in  the  near-bottom  flow  and  the  large  scale  vorticity  variance  in  wave  numbers 

(0.1  <kz  <  1). 

The  velocities  presented  in  Figures  15  - 19  represent  72  second  means  which  were 
sampled  every  other  minute  seventeen  times  (~  33  minutes)  every  other  hour.  The  temperatures 
include  data  from  thermistors  1, 2, 3,4, 5, 6  and  8.  They  are  included  here  to  document  the  absence 
of  stratification.  The  thermistors  present  data  from  0.33  meters  above  bottom  to  2.97  meters  above 
bottom.  The  standard  deviation  in  each  72  second  record  are  offset  and  plotted  with  the  means. 


Vorticity  variance  within  each  72  second  burst,  as  defined  in  Section  5,  is  included  as 
Figures  20  -  24. 
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&  Across  Strait  (dots)  Flow  (m/s)  at  2.4mab  (1 .1  mab  dashed) 


Along  Strait  (solid)  &  Across  Strait  (dots)  Flow  (m/s)  at  2.4mab  (1 .1mab  dashed) 


Figure  16. 


Along  Strait  (solid)  &  Across  Strait  (dots)  Flow  (m/s)  at  2.4mab  (1.1  mab  dashed) 


Figure  17. 


Along  Strait  (solid)  &  Across  Strait  (dots)  Flow  (m/s)  at  2.4mab  (1.1  mab  dashed) 


Figure  19. 


Vorticity  Variance  (X) 


Vorticity  Variance  (X) 


Figure  24, 


7.  CONCLUSIONS 

It  is  apparent  that  the  spectral  density  of  the  along-path  flow  is  filtered  in  a  manner  which 
is  not  yet  well  understood,  albeit  repeatable.  Evaluation  of  vorticity  spectra  must  remain  restricted 
to  regions  below  this  filtering,  until  such  time  the  filtering  is  understood.  When  flow  is  perpendic¬ 
ular  to  the  acoustic  path,  vortex  shedding  appears  to  contaminate  the  spectra  at  wave  numbers 
consistent  with  the  tubing  diameters  of  the  tripod. 
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APPENDIX  A:  Data  Processing  Software 
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seepress^m 


seepress.m 


%  seepress  ~  script  to  convert  pressure  counts  to  pressure  (psia) 

%  for  sparse  sampling  17  sets  *  514  samples  /set 

nn=l; 


%  dint  when  continuous  sampling  is  made 
dint  =  514; 


for  fn=2:12 

eval([load  v_doy_',num2str(fn)]) 
eval([  load  v _j)ress_lnum2str(fh)] ) 

%  find  break  after  32  minutes  of  sampling 
nbreaks  =  find(diff(doy)  >  0.01); 
nbreaks=  [nbreaks  ;length(doy)];  . 
nbegin=l; 


for  nendi=l:length(nbreaks)  %  for  each  32  minute  series 
nend=nbreaks(nendi); 

%  if  complete  32  minute  cycle  compute  3  avgs  for  each  observation 
for  nstart=nbegin:dint:nend 
nstop=nstart+dint-l ; 
if(nstop  >=  nend) 

nbegin  =  nbreaks(nendi)+l; 

bre^; 

end 

if  (nstop  >  length(doy)) 

nstop  =  length(doy); 

disp([last  sample  only  avgs  lnum2str(nstop-nstart+l)]) 

end 

%  converts  4Hz  pressure  counts  to  pressure  (psia)  for  VORT#3 
% 

pc  =  pcounts(nstart:  nstop); 

p=diff(pcO; 

p=[p(l)  p]; 

eiT=find(p  <  0); 

p(err)=p(err)+2'^4; 

frq  =  1000/140; 

N=32  /frq; 

T=(N*le6).  /p; 

T2=T."2; 

T0=30.44344; 

T02=T(r2; 

C=— 4S58.177; 

D=0.066768; 

ddt=l-(T02.  /T2); 

press=C.  *ddt.*(  1-D.  *ddt); 

pstd(nn>=std(press); 

pmean(nn)=mean(press); 

mdoy(nn)=mean(doy(nstart:nstop)); 

nn=nn+l 

end  %  end  of  nstart  loop  (slip  through  hours  dint  at  a  time) 
nbegin=nend+l; 
if  (nbegin  >=  niax(nbreaks)) 
break 

end 

end  %  nbreak  (nendi) 

end  %  end  of  each  v  vel  N  file 
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gettemp.m 


gettemp.m 


%  gettemp  —  to  convert  thermistor  counts  from  bin2mat  to  degrees  C 
% 


load  tempckV3 


%  ck  in  tempckVS: 
%ck  = 

% 


% 

0.0237 

-0,2758 

5.6379 

-0.0182 

% 

0.0238 

-0.2763 

5.6439 

-0.0453 

% 

0.0239 

-0.2770 

5.6352 

-0.0009 

% 

0.0249 

-0.2837 

5.6509 

0.0098 

% 

0.0251 

-0.2852 

5.6561 

-0.0196 

% 

0.0250 

-0.2846 

5.6537 

-0.0081 

% 

0.0252 

-0.2863 

5.6613 

-0.0029 

% 

0.0247 

-0.2822 

5.6451 

-0.0048 

ck2=ck''; 
for  fn=5:15 

eval([load  v__doy_',num2str(fn)]) 
eval([load  v_temp_',num2str(fh)]) 
ntherms=length(temp)  /length(doy) 

%  for  vort3v30  if  (fn  =  4)  %  problem  at  end  of  file 
%  for  vort3v30  ntherms=length(temp)  /(length(doy)+l); 

%  for  vort3v30  end 

nr  =  length(temp)  /ntherms; 
tempv=reshape(temp,ntherms,nr)  VI 0000 ; 
clear  temp 
for  ist=l:8 

temp(:,ist)  =  ck2(14st).*tempv(:,ist)."3+ck2(2,ist).*tempv(:,ist).V+ck2(3,ist).*tempv(:,ist)+ck2(4,ist); 

end 

%if  (fh  =  4) 

%  temp((length(doy  )+!),:)=[]; 

%  end 

eval(["save  v_temp_",num2str(fh)/  temp"]) 

end 


May  15  16:46  1998 
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nxvortb.m 


fixvortb.m 


%  fixvort  —  to  convert  0—2^16  to  —1.5  — >  +1.5  m  /s  ^ 

%  to  read  output  from  bin2mat 

%  jfredericks@whoi.edu  7  122  /94 

nvels=length(velocity)  /length(doy) 

nr=length( velocity)  /nvels 

vel=Teshape(  velocity, nvels,nr)'; 

for  ns=l:nvels 

nn=find(vel(:,ns)  >  32768); 
if  (length(nn)) 

vel(nn,ns)=vel(nn,ns>-65536; 

end 

end 

vel  =  vel.*(1.5  732678);  %  convert  to  vel  (m  Is) 


/ 


May  30  10:20  1995 
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getvelV3.m 


getvelV3.m 


%  procedure  to  convert  VORT#3  velocity  data  to  VORT#2  orientation 
%  upon  completion  real(w„.r)  corresponds  to  X  axis  of  VORT2 
xtop=-vmean(:,3); 

ytop=(vmean(:,ll)-“Vmean(:,7))  /(2*cos(30*pi  /180)); 
xtopcheck=(vmean(:,ll)+vmean(:,7)),  /(2*sin(30*pi  /180)); 

xbot=vmean(:,l); 

ybot=(vmean(:,5>“Vmean(:,9)).  /(2*cos(30*pi  /1 80)); 
xbotcheck=-(vmeaLn(:,5)+vmean(:,9)).  /(2*sin(30*pi  /180)); 

wtop=xtop+i  *y  top; 
wbot=xbot+i*ybot; 

%degrees=45;  %  for  vort3v30.dat 
%degrees=60;  %  for  vort3.617 

%degrees=-56;  %  for  vrt3pu.dat  if  using  3/1  or  0  if  using  9/11 
%rotated 

%topr=wtop.  *thetar; 

%botr=wbot.  *thetar; 


May  30  10:20  1995 
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APPENDIX  B:  Rotation  of  ABC  vorticity  planes  to  XYZ 


Appendix  B.  Rotation  of  Vorticity  Planes  Into  Along-Strait  (x),  Across-Strait(y)  and  Vertical  (z) 
Coordinates 

Assuming  equidistant  length  between  all  sensor  planes  within  each  circulation  path  and  each  sensor 
within  all  horizontal  planes,  each  sensor  can  be  described  as  a  point  source  as  follows: 


(Xb,Yb,Zb)-  (0,0,0) 

(X„Yr,Zr)  =  (L,0,0) 

(Xo,Yo,Zo)  = 

(Xy,Yy,  Zy)  =  ^0,  L, 

(Xw,Yw,Zw)  =  ^L,^L,^L^ 


Each  vorticity  vector  can  be  defined  as  follows: 

a  =  vector  of  unit  length  in  A  direction  =  (ax,ay,az) 
b  =  vector  of  unit  length  in  B  direction  =  (bx,  by,  bz) 
c  =  vector  of  unit  length  in  C  direction  =  (cx,  Cy,  Cz) 

It  can  be  shown  that 


R  = 


1 

2 

1 

2 


-  1 


Ct^xyz 


and 

—  R  *  ^ahc 
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APPENDIX  C:  Data  Logging  &  Unpacking  Software 
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biii2mat.c 


bin2mat.c 


/*  binlmat  -  read  VORT^S  or  VORTil  binary  data  files 

&  output  MATIAB  files  of  time  Ipressure  Ivel  Itemp  *  / 

/*  jjfredericks  Woods  Hole  Oceanographic  Institution 
jfredericks@whoi.edu  6  194  *  / 

/*  NOTE:  must  recompile  with  FLAGS  set  for  V0RT3  or  V0RT2  *  i 

#include  <stdio.h> 

#iiiclude  <strings.h> 
tinclude  <malloah> 

#iiiciude  "math”  /*  symbolic  link  to  /usr  /local  /matlab  /extern  /...  *  / 

tdefine  FILESIZE  229376  /*  No.  bytes  per  disk  dump*  / 

#defuie  MAXFILES  1481  /*  maximum  number  of  disk  dumps  *  / 

#defiDe  ENDBUFF  256  /*  will  not  begin  a  data  record  beyond  FILESIZE  —  ENDBUFF  *  / 

fdeflne  MAXREC  130000  /*  514*17*12^104856  added  extra  */ 

#deflne  RECLENV3  58  /*  VortHS  */ 

#deniie  RECLENV2  74  /*  Vort^l  */ 


#defiiie  FLAGA  Oxab 
#denne  FLAGB  0x53 
/*  ^define  FLAGS  0x52  for  VORTn  *  / 

fdefine  NTEMEWORD  6  /*  number  of  bytes  in  time  word  *  / 

/*  monthly  table  for  computation  of  dayjof_year  *  / 
unsigned  short  monthday[]={0, 31,60,91,121,152, 182,213,244,274,305,335); 

unsigned  char  time[NTIMEWORD],oldsec; 
unsigned  long  nrec; 
int  nvels,ntherms,namex; 

double  doy[MAXREC],*velptr,*temptr,*vels, *temps,pitch[MAXREC],  roll[MAXREC], 
compass[MAXREC],  pcounts [MAXREC]; 
float  fsec,stoptime,startime,atof(),tbisdoy,lastdoy; 

MATFUe  *fp; 

Matrix  *a,*b; 

char  file_out[40]; 


main(argc,argv) 
int  argc; 
char  **argv; 

{ 

unsigned  char  dread  [FILESIZE]; 
int  nb,ic,dsize,nc,status; 


/*  fsec  is  for  the  addition  of  the  140  ms  sampling  rate  between  each  second  *  / 


fprintf(stderr, 

”\n\n\n  TATTLETALE  BINARY  TO  MATLAB  BINARY  CONVERSION  9/l/94\n\n”); 


startime=0; 

stoptime=500.0; 

fprintf(stderr,”sizeof  double:  %ld\n”,sizeof{double)); 

if  (argc  <  2) 

{  fprintf(stdeiT, 


main 
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biii2mat.c 


bin2mat.c 


...main 

”\nUSAGE:  bin2mat  {v2iv3}  [startime  stoptime]  <packed_datafile\n\n'’); 
exit(l); 

} 


if  (strcmp(*++argv,”v2")=0) 

{ 

fprintf(stderr, 

’’Reading  data  in  VORT#2  format/24  acoustic  paths/4  thermistorsXn”); 
dsize=RECLENV2; 
ntherms=4; 
nvels=24; 

1 

else  if  (strcmp(*argv,”v3”)==0) 

{ 

dsize=RECLENV3; 

ntherms=8; 

nvels=12; 

fprintf{stderr, 

"Reading  (%d  byte)  datablocks  in  VORT#3  format/%d  acoustic  paths/%d  thermistorsXn", 
dsize,nvels,ntherms); 

) 

else 

{  fprintf(stderr, 

"XnUSAGE:  bin2mat  {v2lv3}  [startime(doy)  stoptime(doy)]  <packed_datafiie\n\n"); 
exit(l); 

) 

if  (argc  >  2) 

{ 

startime=atof ( *++arg  v); 

fprintf(stderr,"  Start  time  specified:  %g\n", startime); 
if  (argc  >  3) 

{ 

stoptime=  atof(*++argv); 

fprintf(stdeiT,’’  Stop  time  specified:  %g\n",stoptime); 

} 

} 

lastdoy=startime; 

/*  allocate  space  for  temp  &  vel  variable  ncols  *  / 
vels  =  (double  *)  calloc(nvels*MAXREC,sizeof(doable)); 
temps=  (dooble  *)  calloc(ntherms*MAXREC,sizeof(double)); 
velptr  =  vels; 
temptr  =  temps; 

fprintf(stderr,"vels:  %d/  velptr:  %d\n”, vels, velptr); 
fprintf(stderr,"  temps:  %d/  temptr:  %d\n’’,temps,  temptr); 
if  (vels  =  0  II  temps  =  0) 

{ 

fprintf(stderr," trouble  allocating  memory  at  rec:  %d  \n",nrec); 
exit(l); 

} 


nrec  =  0; 
namex  =  0; 

oldsec=  (onsigned  char)  0; 

for  (nb=0;  nb  <  MAXFILES;  nb++) 

{ 

/*  read  in  a  disk  block  *  / 

status  =  fread(dread,sizeof(char),  FILESIZE,  stdin  ); 
if  (status  <=  0) 

{  fprintf(stderr, 

"Error  searching  for  rec  %d  nblock  =  %d  status  =  %d\n". 
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biii2mat.c 


bin2mat.c 


...main 

nrec,nb,status); 

) 

/*  find  beginning  of  valid  records  *  / 
for  (ic=0;  ic  <  FILESIZE;  ic++) 

{ 

if  (((dread[ic]  =  FLAGA)  &&  (dread[ic+l]  =  FLAGB))  && 

((dread[ic+dsize]  ==  FLAGA)  &&  (dread[ic+dsize+l]  =  FLAGB))) 

break; 


if  (ic  >  HLESIZE  -  ENDBUFF)  /*  no  records  in  block  *  / 

{ 

fprintf(stderr,"  LOST  IN  FILE  at  block:  %d\n",nb); 

nrec-=:l; 

dodumpO; 

exit(0); 

) 

if  (status  >=  dsize) 

^rintf(stden‘, 

’’found  first  byte  at  %d  in  block  %d\n'',ic,nb); 
ic++;  /*  found  valid  data  now  skip  over  flag  Oxab  *  / 


/*  BEGIN  getting  RECORDS  from  BLOCK  *  / 

while  (ic  <  FILESIZE-ENDBUFF) 

{ 

if(status  >=  dsize) 

{ 

/*  read  time:  MO  DA  HR  MN  SEC  ff  *  / 
for  (nc=0;  nc  <  NTIMEWORD;  nc++) 
time[  nc]=dread[ ++ic] ; 

if(tiine[4]  =  oldsec)  /*  fix  time  for  140  ms  sampling  rate  *  / 
fsec+=0.140; 

else 

{ 

fsec=0.0; 

oldsec=tinie[4]; 

) 

thisdoy=doy[nrec]=(doublc)  (monthday[time[0]-l]+time[l]+time[2]  /24.0+ 
tiine[3]  /1440.0+((doDble)  time[4]+fsec)  /86400.0); 


ic++;  /*  advance  pointer  to  read  next  byte  *  / 

for  (nc=0;  nc  <  ntherms;  nc++) 

{ 

if(thisdoy  >=  startime) 

*temptr++=(dooble)  (dread[ic++]+dread[ic++]  *256); 

else 

{ic++;ic++;} 

) 

pitch[nrec]=(  (doable)  (dread[ic++]+dread[ic++]*256)-49152)  716384*45; 
roll[nrec]=:(  (doable)  (dread [ic++]+dread[ic++] *256-49152)  716384*45); 

for  (nc=0;  nc  <  nvels;  nc++) 

{ 

If  (thisdoy  >=  startime) 

♦velptr++=  (double)  (dread[ic++]+dread[ic-H-]*256); 
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bin2mat.c 


bin2mat.c 


...main 


else 

{ic++;ic++;} 

) 

compass! nrec]=  (double)  (dread[ic++]+dread[ic-H-]*256); 
pcounts[nrec]=(double)  (dread[ic++] *65536+dread[ic-H-]  *1677721 6+dread[ic++]+dread[ic++] *256); 

/*  force  dump  of  stored  records  if 

reached  stoptime  ...  or  ... 
reached  a  L5  day  gap  ...  or  ... 
read  the  last  fde  block  *  / 

if(lastdoy  >=  stoptime) 

{ 

dodumpO; 

fprintf(stderr,"gone  beyond  stoptimeXn"); 
exit(l); 

} 

if((thisdoy“lastdoy)  >1.5  &&  thisdoy  >  startime) 

dodumpO;  ^  *  / 

if(lastdoy  >=  startime) 
nrec++; 

if(nrec  >=  MAXREC  &&  thisdoy  >  startime) 

{ 

dodumpO; 

exit(O); 

} 

}  /*  end  of  if  (status  >=  dsize)  *  / 

if(status  !=  FILESIZE) 

/*  ran  out  of  file  *  / 

{ 

if(nrec)  dodumpO; 
fprintf(stdeiT, 

’’reached  partial  record  so  quit  readingVn”); 
exit(O); 

} 


fprintfstderrf [dread[%d]:  %0x  ...  %OxWfic,dread[ic],dread[ic-\-l]); 
lastdoy=thisdoy; 

if(ic  >  355  &&  ic  <  365)  /*  hack  a  discovered  problem  *  / 

{ 

fprintf(stderr,” hacked  code  at  nb:%d  \nrec:  %d\n", 
nb,nrec); 
ic  =  ic— 2; 

} 

if(dread[ic++]  1=  FLAG  A  &&  dread[ic]  !=  FLAGB) 

break;  /*  check  for  end  of  data  in  block  *  / 
};  /*  finished  with  data  part  of  disk  block  *  / 


}  /*  finished  with  disk  block  *  / 

fprintf(stderr,”NBLOCKS  READ:  %d,  last  record  contains  %d  bytes\n”,nb,nrec); 
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}  /*  End  of  Main  Loop  *  / 


/♦  routine  to  dump  to  datcfiles  *  / 
dodumpO 

{ 

int  ic; 

if((thisdoy-lastdoy)  <=  1.5)  /*  not  at  multiday  gap  *  / 

nrec  +=  1;  /*  convert  index  to  number  of  records  *  / 

/*  don't  need  to  if  new  day  because  we  don't 
want  to  dump  the  last  record  stored  *  / 
printf(”month(%d)/day(%d)/hr(%d)/min(%d)/sec(%d)  (%d  records)\n'\ 
time[0],tiine[l],time[2],time[3],time[4],nrec); 

a=inxCreateFull(nrec,  1  ,RE  AL); 

sprixltf(file_out,*’%c_doy_%d.mat''/v^namex); 

i^matOpen(fiie_out,’*w” ); 

memcpy(mxGetPr(a),doy,nrec*si2cof(doubIe)); 

mxSetName(a,”doy " ); 

matPutMatrix(fp,a); 

matClose(fp); 

mxFreeMatrix(a); 

sprintf(file__out,’*%c_vel_%d.mat*’/v",namex); 
i^matOpen(file_out,’' w” ); 
a=mxCreateFull((nrec)  *nvels,  1  ,RE  AL); 
mxSetName(a,'' velocity” ); 

memcpy  (inxGetPr(a),  veis,nrec  *n  vels  *sizeof(dou  ble)) ; 

matPutMatrix(ip,a); 

matClose(^); 

mxFreeM  atrix(a); 

sprintf(file_out,”%c_temp_%d,mat”,  'v  ",namex); 
f^matOpen(file_out,”  w” ); 
a=inxCreateFull((nrec)  *ntherms,  1  ,RE  AL); 
mxSetName(a,"temp” ); 

memcpy  (mxGetPr(a),  temps, nrec  *n  therms  *slzeof(doubIc)); 

matPutMatrix(ip,a); 

matClose(fp); 

mxFreeM  atrix(a) ; 

sprintf(file_out,"%c__press_%d.mat",  'v  ",namex); 

^matOpen(file_out, "  w" ) ; 
a=mxCreateFull(nrec,  1  ,RE  AL); 
inxSetName(a,”pcounts" ); 

memcpy(mxGetPr(a),pcounts,nrec*sizeof(double)); 

matPutMatrix(fp,a); 

matClose(fp); 

mxFreeMatrix(a); 

sprintf(file_out,”  %c_prc_%d.mat”,  V  ',namex); 
fp=matOpen(file_out,  ”  w" ) ; 
a=mxCreateFull(nrec,  1  ,RE  AL); 
mxSetName(a,  "pitch” ); 

memcpy(mxGetPr(a),pitch,nrec*sizeof(double)); 

matPutMatrix(fp,a); 

mxSetName(a,"roir' ); 

memcpy(mxGetPr(a),roll,nrec*si2eof(doobIe)); 

matPutMatrix(fp,a); 

inxSetName(a,  "compass" ); 

memcpy(mxGetPr(a),compass,nrec*si2eof(double)); 


dodump 
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...dodump 

matPutMatrix(fp,a); 

matClose(fp); 

mxFreeMatrix(a); 

namex++; 

if((thisdoy-lastdoy)  <=  1.5)  /*  not  a  inter  day  gap  *  / 

{ 

lprintf(stderr/’nrec:  %d\nvelptr:  %ld\ntemptr:  %ld\n", 
nrec,velptr,temptr); 
nrec=^-l; 
velptr  =  vels; 
temptr  =  temps; 

) 

else  /*  copy  last  to  first  for  next  day  *  / 

{  . 

doy[0]=:doy[nrec]; 
pitch[0]=pitch[nrec] ; 
roll[0]=roll[nrec]; 
compass[0]=compass[nrec] ; 
pcounts[0]  =pcounts  [nrec] ; 

for  (ic=0;  ic  <  nvels;  ic++) 
velptr — ; 

for  (ic=0;  ic  <  ntherms;  ic++) 
temptr — ; 

memcpy  (vels,  velptr, nvels  *sizeof(double)); 
memcpy(temps,temptr,ntherms*sizeof(doiible)); 

fprintf(stderr,"nrec:  %d\n velptr:  %ld\ntemptr:  %ld\n", 
nrec, velptr, temptr); 


nrec=0; 
velptr  =  vels; 
temptr  =  temps; 
for  (ic=0;  ic  <  nvels;  ic++) 
velptr++; 

for  (ic=0;  ic  <  ntherms;  ic++) 
temptr++; 

fprintf(stderr,’'nrec:  %d\nvelptr:  %ld\ntemptr:  %ld\n", 
nrec,  velptr, temptr); 

) 

}  /*  end  of  dump  to  matlab  files  *  / 
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10  PRINT  "  V0RT3  BASS,  26-JUNE-94  ROM  VERSION  from  IMCSBASS  17-March-94" 

11  PRINT  "  514  RECORDS  EVERY  2  MIN  FOR  16  BLOCKS  (32  MIN)  EVERY  2  HR," 

12  PRINT  "  EVERY  3  DAYS,  EVERY  OTHER  MONTH." 

'20  Q=&H4000:REM  DATAFILE  STARTS  AT  4000  IN  RAM 

20  Q=&H0143:REM  DATAFILE  STARTS  AT  0143  IN  ROM 


'F=PWR  FLAG,  0  FOR  NOTHING  ON,  1  FOR  P&S,  2  FOR  P,S,&0. 

'D=DO  CALIBRATE  A/D  AND  READ  COMPASS  LATCH 

'ACMs,OBSs,PRES,&  SEABIRDS  MEASURED  FOR  15  MINUTES  EVERY  HOUR  AT  4  Hz. 
'TO  LATCH  PWR,  FIRST  PWR  ON  AND  WASTE  TIME  WITH: 

'PSET  16: CALL  &HFFD9 ,0:CALL  &HFFD9,0,  THEN: 


'COMPASS  ON  PSET  15 

'PRESS  ON  PSET  7,14:PCLR  7,16,14 

'SEABIRDS  ON  PSET  6,14:PCLR  6,16,14 

'OBSs  ON  PSET  2,14,:PCLR  2,16,14. 

'TO  PREVENT  TURNING  THESE  OFF  WHEN  SELECTING  AUX  WITH  D\6, 

'HOLD  D\4  HI  WITH  PSET  16: CALL  &HFFD9 , 0 : CALL  &HFFD9,0 

'PSET  1,14:PCLR  1,16  WHERE  CALL  &HFFD9  IS  A  2.2MS  TEST  OF  UART  BUFFER. 
100  ASM  &H9E,DB  &H00:REM  CONTROL  DISABLE 

102  SLEEP  0:PCLR  0, 1 , 2 , 3 , 4 ,5 , 6 , 7 ,8 , 14 , 15 :PSET  14,16:REM  ESTABLISH  DDRs 
108  SLEEP  240:PCLR  14, 16: REM  LET  A/D  CALIBRATE  FINISH  IF  STARTED. 

110  PRINT  "  TYPE  WAKE  TO  GET  CONTROL  (5sec)" 

112  X  =  0:  STORE  X,#4,&H01010101  :  X  =  0:  ITEXT  X,500 
114  X  =  0:  IF  GET(X,#4)<>&H57414B45  GOTO  300 


130  RTIME 

131  PRINT  #02, "OLD  TIME=" ,?( 3 ),"/",?( 4 ),"/",?( 5 ), "  ",?( 2 1 0) 

132  INPUT  "  SET  TIME?  1=YES",J 

133  IF  JO  1  GOTO  160 

134  INPUT  "  THE  YEAR  0-99  ",?(5) 

136  INPUT  "  THE  MONTH  1-12  ",?(4) 

138  INPUT  "  THE  DAY  1-31  ",?(3) 

140  INPUT  "  THE  HOUR  0-23  ",?(2) 

142  INPUT  "  THE  MINUTE  0-59  ",?(!) 

144  INPUT  "  THE  SECOND  0-59  ",?(0) 

150  STIME 

160  RTIME 

162  PRINT  #02,"NEW  TIME=",?(3),"/",?(4),"/",?(5),"  ",?( 2) ?(1 0) 
164  INPUT  "  STOP?  1=YES",J 
166  IF  J<>  1  GOTO  300 
230  STOP 


300  ONERR  100 
1000  DIM(1024,8192) 

'  ASSEMBLY  ROUTINES 

1010  X=0 : A=0 : B=0 : C=0 : D=0 : E=514 : F=0 : G=0 : H=0 : 1=0 : J=0 : K=0 : L=0 : M=0 : N=0 : 0=0 

1015  PRINT  "  A  B  C  D  E  F  G  H"; 

1016  PRINT  "  I  J  K  L  M  N  O  P  T  X" 

1020  GOSUB  8000 ;REM  FIRST  PASS 

1025  PRINT  #5H,A,B,C,D,E,F,G,H,I,J,K,L,M,N,0,P,T,X 

1030  GOSUB  8000 ;REM  SECOND  PASS 

1031  PRINT  #5H,A,B,C,D,E,F,G,H,I,J,K,L,M,N,0,P,T,X 
1050  N=0:REM  INITIALIZE  COUNTER 

1080  GOTO  2000 


'SUBROUTINE  TO  TURN  ON  PRESSURE  AND  SEABIRDS 

1200  PSET  16:CALL  &HFFD9 , 0 : CALL  &HFFD9,0:REM  LET  PWR  SETTLE 
1210  PSET  6,7,14:PCLR  6 , 7 , 16 , 14 : F=1 :RETURN 
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'SUBROUTINE  TO  TURN  ON  OBSS,  PRESSURE,  AND  SEABIRDS 

1230  PSET  16: CALL  &HFFD9 , 0 : CALL  &HFFD9,0;REM  LET  PWR  SETTLE 

1240  PSET  2,6,7,14:PCLR  2 , 6 , 7 , 16 , 14 :F=2 :RETURN 

'SUBROUTINE  TO  TURN  OFF  EVERYTHING 

1250  PCLR  2, 6, 7, 14: PSET  14:PCLR  14 : F=0 : RETURN 

'SUBROUTINE  TO  ENABLE  AUXILLIARIES  ON  MULTIPLEXOR 
1260  PSET  16: REM  LET  PP9R  SETTLE 

1270  PSET  1,14:PCLR  1 , 1 6 : RETURN : REM  DISABLES  POWER  LATCH 
'  CALIBRATE  A/D 

2000  SLEEP  0:PSET  15: SLEEP  600: PSET  16: REM  PWR  ON, CLEAR  LATCHED  PWR 

2050  CALL  T,0:REM  FIND  COMPASS  GAP 

2060  Z=COUNT( 40 )-l: SLEEP  50:REM  Z=aZiinuth, COMPASS 

'2065  PRINT  "GOT  COMPASS  GAP" 

2070  PCLR  15:SLEEP  240:REM  A/D  CALIBRATE; 2, 882, 040  CYCLES  AT  1.2288  MHz 

2080  PCLR  14, 16: REM  POWER  OFF  AND  PREPARE  TO  LATCH  PWR 

2082  D=0:REM  CLEAR  LATCH  TO  DO  A/D  AND  Zn 

2090  GOSUB  1200: REM  TURN  ON  SEABIRD , PRESS 

2100  SLEEP  0:P=SDI(24) ;REM  SAMPLE  AND  READ  PRESSURE 

3000  RTIME:IF  E>0  GOTO  3040 

3001  IF  ?(4)%2=0  GOTO  3030:REM  RUN  IN  JA,MAR,MAY, JUL,S  ( JUL,S,N, JA,MAR) 
3006  IF  ?(3)%3<>1  GOTO  3030 :REM  RUN  EVERY  THIRD  DAY 

3008  IF  ?(2)%2>0  GOTO  3030 :REM  ONLY  RUN  EVERY  OTHER  HOUR 
3010  IF  ?(1)>32  D=l:GOTO  3030:REM  OFF  MINUTE  33  TO  END  OF  HOUR 
3012  IF  ?(0)>0  GOTO  3030 :REM  CONTINUE 

3014  IF  ?(1)%2=0  E=514 

3015  IF  D=1  GOTO  2000 ;REM  COMPASS  AND  A/D  ONCE/HR 
3020  GOTO  2090: REM  POWER  ON  PRESSURE  AND  START  CYCLE 

3030  GOSUB  1250: SLEEP  100: GOTO  3000 :REM  EVERYTHING  OFF 

3040  E=E-1:X=0:REM  INITIALIZE  DATAFILE 

3050  STORE  X, #2 , &HAB53 :REM  KEY  WORD  AB  V3  FOR  VORT3  OR  AB53 

3051  STORE  X,#l,?(4) :REM  MONTH 

3052  STORE  X, #1 , ? ( 3 ) :REM  DAYS 
3055  STORE  X, #1 , ? (2 ) :REM  HOURS 
3060  STORE  X,#l,?(l) :REM  MINUTES 
3070  STORE  X,#1,?(0) :REM  SECONDS 
3080  STORE  X,#1,N:REM  COUNTER 

3100  C=&H73000000+Q+X:REM  MUX  ADDRESS  AND  DATAFILE  POINTER 

3110  CALL  A,C,X:REM  CALL  A/D  ROUTINE, ACM  THEN  AUX, POWER  OFF  AT  END 

3111  X=X%&H10000-Q:REM  GET  BACK  THE  DATAFILE  POINTER 

'  SUBTRACT  ROUTINE  (CALL  B)  EXPECTS  FLAGGED  DATA  AT  6F38. 

3119  C=&H73000000+Q+X:REM  DATAFILE  LOCATION 

3120  CALL  B,C,X:REM  CALL  SUBTRACT  AND  TRANSFER 

3121  X=X%&H10000-Q:REM  GET  BACK  THE  BASIC  DATAFILE  POINTER 
3125  STORE  X,#2,Z:REM  COMPASS 

3130  STORE  X,#4,P:REM  PRESSURE 

'3170  X=0:REM  REINITIALIZE  DATAFILE 
'3180  PRINT  #4H,GET(X,#2) :REM  KEYWORD 

'3185  PRINT  #02,GET(X,#1),"  " ,GET(X, #1 ) , "  ";:REM  MONTH  DAY 

'3190  PRINT  #02,GET(X,#1) ,":",GET(X,#1) ,":",GET(X,#1) ,"  " ; :REM  HR:MM:SS 

'3200  PRINT  #02H,  GET(X,#1) ;REM  COUNT 

'3205  FOR  L=1  TO  10:REM  T(0:7) , PITCH, ROLL 


58 


'3206  PRINT  #04H,GET(X,#2)  ; 

'3207  NEXT  L: PRINT 

'3210  FOR  M=1  TO  3:REM  1  VORT  SENSOR 
'3220  FOR  L=1  TO  4: REM  4  WORDS/LINE 

'3230  PRINT  #04H,GET(X,#2) ";;REM  EACH  SENSOR  ON  A  LINE 
'3240  NEXT  L: PRINT: REM  FOUR  AXES/LINE 
'3250  NEXT  M 

'3262  PRINT  #3  ,GET(X,  #2)  ,  "  REM  COMPASS 
'3264  PRINT  #08H,GET(X, #4 ) , "  ";;REM  PRESSURE 
'3275  PRINT: REM  SPACE 
3280  PRINT  {0,57};:REM  TAKES  59  ms. 


3290  N=N+1 

3295  SLEEP  14 : P=SDI ( 24 ) :REM  SAMPLE  COUNTERS  AND  READ  PRESSURE 
3300  GOTO  3000 


8000 

X=&H7300: 

REM  MULTIPLEXOR  LIST 

8010 

ASM 

X,DW 

&H0000;DW 

&H0101;DW 

&H0202;DW 

&H0303:REM 

ACM 

1/T1-T4 

8015 

ASM 

X,DW 

&H0808;DW 

&H0909;DW 

&H0A0A;DW 

&H0B0B:REM 

ACM 

1/T5-T8 

8020 

ASM 

X,DW 

&H208A;DW 

&H218B;DW 

&H22FF;DW 

&H23FF:REM 

ACM 

2/PITCH 

8025 

ASM 

X,DW 

&H28FF;DW 

&H29FF;DW 

&H2AFF;DW 

&H2BFF 

8030 

ASM 

X,DW 

&H4000;DW 

&H4100;DW 

&H4200;DW 

&H4300:REM 

ACM 

3 

8035 

ASM 

X,DW 

&H4800;DW 

&H4900;DW 

&H4AOO;DW 

&H4B00 

'8040  ASM  X,DW  &H6000;DW  &H6100;DW  &H6200;DW  &H6300:REM  ACM  4 

'8045  ASM  X,DW  &H6800;DW  &H6900;DW  &H6A00;DW  &H6B00 

'8050  ASM  X,DW  &H8000;DW  &H8100;DW  &H8200;DW  &H8300:REM  ACM  5 

'8060  ASM  X,DW  &H8800;DW  &H8900;DW  &H8A00;DW  &H8B00 

'8070  ASM  X,DW  &HA000;DW  &HA100;DW  &HA200;DW  &HA300:REM  ACM  6 

'E8080  ASM  X,DW  &HA800;DW  &HA900;DW  &HAA00;DW  &HAB00 

8090  ASM  X,DW  &HFF00;DW  &HFF00:REM  END  OF  LIST 
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8100  X=&H3800:A=X:REM  A/D  ROUTINE 

'  X  IS  POINTER  TO  AUX  ARRAY,  AB  IS  POINTER  TO  FLAGGED  CONVERSION  STORAGE 
8105  ASM  X,STD  &H7000:REM  POINTER  TO  DATAFILE 

8107  ASM  X,LDD  #&H6F38;STD  &H7002;REM  POINTER  TO  FLAGGED  CONVERSIONS 
8110  ASM  X,OIM  &H04,&H17;PSHX;LDX  #&H0600:REM  POWER  ON,  WAIT  5ms 
8120  D=X:ASM  X,DEX;BNE  D;PULX:REM  WAIT  (4*.81us/LOOP) 

8150  ASM  X,CLR  &H11;LDAA  #&H10;STAA  &H10:REM  DISABLE  UART,  SET  SCI  FOR 
8160  ASM  X,OIM  8,&H11:REM  INTERNALLY  CLOCKED  8  BIT  DATA  RECEIVE  ENABLE 
'DUMMY  PULSE  SENT  OUT  TO  INITIALIZE 

8170  ASM  X,OIM  &HAB,&H15;AIM  &HBF,&H15:REM  DUMMY  AUX  LOCATION 

8171  ASM  X,OIM  &H80,&H03;OIM  &H80,&H03;OIM  &H80,&H03 ;OIM  &H80,&H03 
'START  TIMING  P27=l 

8172  ASM  X,AIM  &H7F, &H03 :REM  REMOVE  PULSE  P27=0 

8174  D=X:ASM  X,LDAB  #&H02;BITB  &H03;BEQ  D:REM  CHECK  P21  FOR  DONE 

8179  C=X:ASM  X,AIM  &H14 , &H15 :REM  ENTRY  POINT  FOR  AUX, MASK  PORT  5 

8180  ASM  X,LDAA  &H00,X;INX;INX:REM  TOP  OF  MUX  LIST, MOVE  TO  NEXT 
8185  ASM  X,ORAA  &H15;STAA  &H15:REM  PUT  MUX  ON  PORT  5 

8200  ASM  X,OIM  &H80,&H0370IM  &H80 , &H03 ;OIM  &H80 , &H03 ;OIM  &H80,&H03 
'START  TIMING  P27=l 

8210  ASM  X,AIM  &H7F, &H03 :REM  REMOVE  PULSE  P27=0 

8220  D=X:ASM  X,LDAB  #&H02;BITB  &H03;BEQ  D:REM  CHECK  P21  FOR  DONE 

8240  ASM  X,LDAB  &H15;ANDB  #&H14:REM  MASKED  PORT  5  IN  B 

8250  ASM  X,LDAA  &H00,X;INX;INX;REM  GET  NEXT  MUX  WORD, ADVANCE 

8255  ASM  X,ABA?STAA  &H15:REM  ADD  MUX  WORD  TO  MASKED  PORT  5,  STORE  P5 

8280  ASM  X,OIM  &H80 , &H03 ;OIM  &H80 , &H03 ;OIM  &H80 , &H03 ;OIM  &H80,&H03 

/ 
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START  TIMING  P27=l 

8290  ASM  X,AIM  &H7F, &H03 :REM  REMOVE  PULSE  P27=0 

8400  ASM  X,CPX  0;CPX  0;CPX  0;CPX  0 

8402  ASM  X,CPX  0;CPX  0;CPX  0;CPX  0 

8403  ASM  X,CPX  0;CPX  0;CPX  0;CPX  0 

8404  ASM  X,CPX  0;CPX  0;CPX  0;CPX  0 

8410  ASM  X,LDAA  «eHll;LDAB  &H12:REM  TRIGGER  A  READ 
8420  ASM  X,CPX  0;CPX  0;CPX  0;CPX  0;REM  16ECLOCK 

8430  ASM  X,LDAA  &H11;LDAB  &H12:REM  HIGH  BYTE  IN  B  AND  READ  NEXT 
8440  ASM  X,CPX  0;CPX  0;CPX  0;CPX  0:REM  16ECLOCK 

8450  ASM  X,LDAA  &H12;REM  LOW  BYTE  IN  A  DON'T  READ  NEXT 
'REM  FLIP  THE  BITS  OF  THE  WORDS 

8451  ASM  X,ROLA;RORB;ROLA;RORB;ROLA;RORB;ROLA;RORB 

8452  ASM  X,ROLA;RORB;ROLA;RORB;ROLA;RORB;ROLA;RORB;ROLA 

8480  ASM  X,TIM  &H20 , &H03 ;BEQ  G:REM  TEST  BOTH  RECEIVED 
8488  ASM  X,XGDX;RORB;BCS  M:REM  BRANCH  IF  AUXILLIARY  (ODD  X) 

8490  ASM  X,ROLB;XGDX;CLRA;CLRB;BRA  G:REM  IF  NOT  RECEIVED , CLEAR 
8495  M=X:ASM  X,ROLB;XGDX 

'  PUT  AWAY  TO  §(50)  ARRAY  THE  FLAGGED  VALUES 

8510  G=X:ASM  X,PSHX;LDX  &H7002;STD  0,X;INX;INX;STX  &H7002;PULX 

8520  ASM  X,LDAB  #&HFF;EORB  &H00,X;BEQ  L;JMP  D:REM  CHECK  FOR  END  OF  LIST 

8550  L=X:ASM  X,LDAB  #&H02:REM  TEST  FOR  HOLD  PULSE 
8555  D=X:ASM  X,BITB  &H03;BEQ  DtREM  CHECK  P21 

8572  ASM  X,CPX  0;CPX  0;CPX  0;CPX  0 

8573  ASM  X,CPX  0;CPX  0;CPX  0;CPX  0 

8574  ASM  X,CPX  0;CPX  0;CPX  0;CPX  0 

8575  ASM  X,CPX  0;CPX  0;CPX  0;CPX  0 

8580  ASM  X,LDAA  &H11;  LDAB  &H12:REM  TRIGGER  A  READ 
8590  ASM  X,CPX  0;CPX  0;CPX  0;CPX  0:REM  16ECLOCK 

8600  ASM  X,LDAA  &H11;LDAB  &H12:REM  HIGH  BYTE  IN  B  AND  READ  NEXT 
8610  ASM  X,CPX  0;CPX  0;CPX  0;CPX  0:REM  16ECLOCK 

8620  ASM  X,LDAA  &H12:REM  LOW  BYTE  IN  A  DON'T  READ  NEXT 

'REM  FLIP  THE  BITS  OF  THE  WORDS  PUTS  HIGH  BYTE  IN  A,  LOW  BYTE  IN  B 

8621  ASM  X,ROLA;RORB;ROLA;RORB;ROLA;RORB;ROLA;RORB 

8622  ASM  X,ROLA;RORB;ROLA;RORB;ROLA;RORB;ROLA;RORB;ROLA 

8630  ASM  X,TIM  &H20 , &H03 ;BEQ  H:REM  TEST  BOTH  RECEIVED 
8640  ASM  X,XGDX;RORB;BCS  N:REM  BRANCH  IF  AUXILLIARY  (ODD  X) 

8645  ASM  X,ROLB;XGDX;CLRA; CURB; BRA  HrREM  IF  NOT  RECEIVED , CLEAR 
8650  N=X:ASM  X,ROLB;XGDX 

'  FINISH  PUT  AWAY  TO  §(50)  ARRAY  THE  FLAGGED  VALUES 

8660  H=X:ASM  X,PSHX;LDX  &H7002;STD  0 ,X; INX; INX;STX  &H7002;PULX 

8670  ASM  X,XGDX;RORB;BCS  0:REM  BRANCH  IF  AUXILLIARY  (ODD  X) 

8680  ASM  X,ROLB;XGDX;OIM  &H40,&H15;OIM  &H01,&H17:REM  PREPARE  FOR  AUX 
8690  ASM  X,LDX  #&H7301;LDD  &H7000;STD  &H7002;JMP  C:REM  DO  AUX 

'  EXIT  AFTER  DOING  ACM  AND  AUX 

8700  0=X:ASM  X,ROLB;XGDX 

8705  ASM  X,CLR  &H11;REM  DISABLE  SCI 

8710  ASM  X,LDAA  #5;STAA  &H10:REM  BAUD  RATE  FROM  TIMERl  (#5  GIVES  9600) 
8720  ASM  X,OIM  &H42 , &H1B ;REM  ENABLE  TIMER2  INTS  AND  E/128 
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'8730  ASM  X,LDAA  #95;STAA  &H1C;REM  TIMER  COUNT  TO  GET  100  Hz 
8740  ASM  X,OIM  &H02 , &H11 :REM  02  RESTART  UART  NO  RECEIVE 
'  reconf  uart,  ENABLE  IT,  power  off  (PHYSICALLY  ENABLE  IT) 
8750  ASM  X,AIM  &H14 , &H15 :REM  PUT  MULTIPLEXORS  ON  PARK 
8760  ASM  X,AIM  &HFA, &H17 :REM  P60,P62=0,AUX  &  POWER  OFF 
8775  ASM  X,LDD  &H7002 
8780  ASM  X,RTS 


9100  B=X:REM  SUBTRACT  AND  TRANSFER  SUBROUTINE 

'  POINTERS;  m(7000)  =  DATA  ARRAY  PERMANENT  STORAGE.  PASSED  IN  AB 
'  m(7002)  =  FLAGGED  AXIS  PAIRS  STORAGE 

'  6F38  IS  §(50)  IN  RAM  MUST  BE  SAME  AS  A/D  ROUTINE 

'  m(  X  )  =  AUX  LIST  FOR  END  OF  LIST  CHECKING.  PASSED  IN  X 


9110 

ASM 

X 

9115 

ASM 

X 

9120 

K=X: 

9140 

ASM 

X 

9150 

ASM 

X 

9152 

ASM 

X 

9160 

ASM 

X 

9170 

ASM 

X 

9200 

X 

II 

Q 

9203 

REM 

i 

9205 

ASM 

X 

:  ASM  X,PSHX;  LDX  &H7002;  LDD  0,X;  INX;INX;REM  GET  WORD 


:  ASM  X,  INX;INX;STX  &H7002 ;PULX:REM  MOVE  ON  TO  NEXT  PAIR 


'  CHECK  FOR  END  OF  LIST.  AUX  LIST=  FOUR  BYTES  /  AXIS 
9240  ASM  X,INX;INX;INX;INX 

9250  ASM  X,LDAB  #&HFF;REM  CHECK  FOR  END  OF  LIST 

9255  ASM  X,EORB  &H00,X 

9260  ASM  X,BNE  K;REM  LOOP  IF  NOT  DONE 

9265  ASM  X,LDD  &H7000;REM  LOAD  DATAFILE  POINTER  FOR  RETURN 

9270  ASM  X,RTS:REM  EXIT 

9275  F=X 
'REM  FLAG  MISSED 
9280  ASM  X,LDD  #&H8000 
9290  ASM  X,BRA  D 

'  FIND  COMPASS  N+1  TRANSMIT  GAP 

9500  T=X:ASM  X,LDAA  &H03:REM  PORT  2 

9510  D=X;ASM  X,CMPA  &H03;BEQ  D;REM  FIND  EDGE 

9520  ASM  X,LDAA  &H03;LDX  #&H007B:REM  Ims  GAP 

9530  P=X;ASM  X,CMPA  &H03;BNE  D;DEX;BNE  P:REM  FIND  Ims  GAP 

9540  ASM  X,RTS 

9900  RETURN 


9999  END 
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10  /JR-INT  ”  VORT2  BASS,  26- JUNE-94  ROM  VERSION  from  IMCSBASS  17-March-94'' 

11  PRINT  "  514  RECORDS  EVERY  2  MIN  FOR  16  BLOCKS  (32  MIN)  EVERY  2  HR," 

12  PRINT  "  EVERY  3  DAYS,  EVERY  OTHER  MONTH." 

'20  Q=&H4000:REM  DATAFILE  STARTS  AT  4000  IN  RAM 

20  Q=&H0143:REM  DATAFILE  STARTS  AT  0143  IN  ROM 


'F=PWR  FLAG,  0  FOR  NOTHING  ON,  1  FOR  P&S,  2  FOR  P,S,&0. 

'D=DO  CALIBRATE  A/D  AND  READ  COMPASS  LATCH 

'ACMs,OBSs,PRES,&  SEABIRDs  MEASURED  FOR  15  MINUTES  EVERY  HOUR  AT  4  Hz . 
'TO  LATCH  PWR,  FIRST  PWR  ON  AND  WASTE  TIME  WITH: 

'PSET  16: CALL  &HFFD9,0:CALL  &HFFD9,0,  THEN: 

'COMPASS  ON  PSET  15 

'PRESS  ON  PSET  7 , 14 :PCLR.  7 , 16 , 14 

'SEABIRDS  ON  PSET  6,14:PCLR  6,16,14 

'OBSS  ON  PSET  2,14,:PCLR  2,16,14. 

'TO  PREVENT  TURNING  THESE  OFF  WHEN  SELECTING  AUX  WITH  D\6, 

'HOLD  D\4  HI  WITH  PSET  16: CALL  &HFFD9 , 0 : CALL  &HFFD9,0 

'PSET  1,14:PCLR  1,16  WHERE  CALL  &HFFD9  IS  A  2.2MS  TEST  OF  UART  BUFFER. 
100  ASM  &H9E,DB  &H00:REM  CONTROL  DISABLE 

102  SLEEP  0:PCLR  0 , 1 , 2 , 3 ,4 ,5,6 ,7 ,8 ,14 , 15:PSET  14,16:REM  ESTABLISH  DDRs 
108  SLEEP  240:PCLR  14,16:REM  LET  A/D  CALIBRATE  FINISH  IF  STARTED. 

110  PRINT  "  TYPE  WAKE  TO  GET  CONTROL  (5sec)" 

112  X  =  0:  STORE  X, #4 , &H01010101  :  X  =  0:  ITEXT  X,500 
114  X  =  0:  IF  GET(X-,#4)<>&H57414B45  GOTO  300 


130  RTIME 

131  PRINT  #02, "OLD  TIME=" ,?( 3 ),"/",?( 4 ),"/",?( 5) , " 

132  INPUT  "  SET  TIME?  1=YES",J 

133  IF  JO  1  GOTO  160 

134  INPUT  "  THE  YEAR  0-99  ",?(5) 

136  INPUT  "  THE  MONTH  1-12  ",?(4) 

138  INPUT  "  THE  DAY  1-31  ",?(3) 

140  INPUT  "  THE  HOUR  0-23  ",?(2) 

142  INPUT  "  THE  MINUTE  0-59  ",?(!) 

144  INPUT  "  THE  SECOND  0-59  ",?(0) 

150  STIME 

160  RTIME 

162  PRINT  #02, "NEW  TIME=" ,?( 3 ),"/",?( 4 ),"/",?( 5) , " 
164  INPUT  "  STOP?  1=YES",J 
166  IF  JO  1  GOTO  300 
230  STOP  < 


",?(2) ,":",?(!) ,":",?(0) 


",?(2),":",?(1),":",?(0) 


300  ONERR  100 
1000  DIM(1024,8192) 

'  ASSEMBLY  ROUTINES 

1010  X=0 :  A=0 :  B=0 :  C==0 :  D=0 :  E=514 :  F=0 :  G=0 :  H=0 :  1=0 :  J=0 :  K=0 :  L=0 :  M=0 :  N=0 :  0=0 

1015  PRINT  "  A  B  C  D  E  F  G  H" ; 

1016  PRINT  "  I  J  K  L  M  N  O  P.  T  X" 

1020  GOSUB  8000 :REM  FIRST  PASS 

1025  PRINT  #5H,A,B,C,D,E,F,G,H,I,J,K,L,M,N,0,P,T,X 

1030  GOSUB  8000 :REM  SECOND  PASS 

1031  PRINT  #5H,A,B,C,D,E,F,G,H,I,J,K,L,M,N,0,P,T,X 
1050  N=0:REM  INITIALIZE  COUNTER 

1080  GOTO  2000 


'SUBROUTINE  TO  TURN  ON  PRESSURE  AND  SEABIRDS 

1200  PSET  16: CALL  &HFFD9 , 0 : CALL  &HFFD9,0:REM  LET  PWR  SETTLE 
1210  PSET  6,7,14:PCLR  6 , 7 , 16 , 14 :F=1 :RETURN 


63 


'SUBROUTINE  TO  TURN  ON  QBSs ,  PRESSURE,  AND  SEABIRDS 

1230  PSET  16: CALL  &HFFD9 , 0 : CALL  &HFFD9,0:REM  LET  PWR  SETTLE 

1240  PSET  2,6,7,14:PCLR  2 , 6 , 7 , 16 , 14 : F=2 : RETURN 

'SUBROUTINE  TO  TURN  OFF  EVERYTHING 

1250  PCLR  2, 6, 7, 14: PSET  14:PCLR  14 :F=0: RETURN 

'SUBROUTINE  TO  ENABLE  AUXILLIARIES  ON  MULTIPLEXOR 
1260  PSET  16:  REM  LET  PWR  SETTLE 

1270  PSET  1,14:PCLR  1 , 16 : RETURN: REM  DISABLES  POWER  LATCH 
'  CALIBRATE  A/D 

2000  SLEEP  0;PSET  15: SLEEP  600: PSET  16: REM  PWR  ON, CLEAR  LATCHED  PWR 

2050  CALL  T,0:REM  FIND  COMPASS  GAP 

2060  Z=COUNT( 40 )-l: SLEEP  50:REH  Z=aZiinuth, COMPASS 

'2065  PRINT  "GOT  COMPASS  GAP" 

2070  PCLR  15: SLEEP  240;REM  A/D  CALIBRATE, *2, 882, 040  CYCLES  AT  1.2288  MHz 

2080  PCLR  14, 16: REM  POWER  OFF  AND  PREPARE  TO  LATCH  PWR 

2082  D=0:REM  CLEAR  LATCH  TO  DO  A/D  AND  Zn 

2090  GOSUB  1200:REM  TURN  ON  SEABIRD , PRESS 

2100  SLEEP  0:P=SDI(24)  :REM  SAMPLE  AND  READ  PRESSURE 

3000  RTIME:IF  E>0  GOTO  3040 

3001  IF  ?(4)%2=0  GOTO  3030:REM  RUN  IN  JA,MAR,MAY, JUL,S  ( JDL,S,N, JA,MAR) 
3006  IF  ?(3)%3<>1  GOTO  3030 :REM  RUN  EVERY  THIRD  DAY 

3008  IF  ?(2)%2>0  GOTO  3030:REM  ONLY  RUN  EVERY  OTHER  HOUR 
3010  IF  ?(1)>32  D=l:GOTO  3030:REM  OFF  MINUTE  33  TO  END  OF  HOUR 
3012  IF  ?(0)>0  GOTO  3030 ;REM  CONTINUE 

3014  IF  ?(1)%2=0  E=514 

3015  IF  D=1  GOTO  2000:REM  COMPASS  AND  A/D  ONCE/HR 
3020  GOTO  2090 :REM  POWER  ON  PRESSURE  AND  START  CYCLE 

3030  GOSUB  1250:SLEEP  100:GOTO  3000:REM  EVERYTHING  OFF 

3040  E=E-1:X=0:REM  INITIALIZE  DATAFILE 

3050  STORE  X, #2 , &HAB52 :REM  KEY  WORD  AB  V2  FOR  VORT2  OR  AB52 

3051  STORE  X,#i',?(4):REM  MONTH 

3052  STORE  X,#ll,?(3)  :REM  DAYS 
3055  STORE  X,#i,?(2) :RBM  HOURS 
3060  STORE  X,#l,?(l)  :REM  MINUTES 
3070  STORE  X,#1,?{0) :REM  SECONDS 
3080  STORE  X,#1>N:REM  COUNTER 

3100  C=&H73000000+Q+X:REM  MUX  ADDRESS  AND  DATAFILE  POINTER 

3110  CALL  A,C,X:REM  CALL  A/D  ROUTINE, ACM  THEN  AUX, POWER  OFF  AT  END 

3111  X=X%&H10000-Q:REM  GET  BACK  THE  DATAFILE  POINTER 

'  SUBTRACT  ROUTINE  (CALL  B)  EXPECTS  FLAGGED  DATA  AT  6F38. 

3119  C=&H73000000+Q+X:REM  DATAFILE  LOCATION 

3120  CALL  B,C,X:REM  CALL  SUBTRACT  AND  TRANSFER 

3121  X=X%&H10000-Q:REM  GET  BACK  THE  BASIC  DATAFILE  POINTER 
3125  STORE  X,#2,Z:REM  COMPASS 

3130  STORE  X,#4,P:REM  PRESSURE 

'3170  X=0:REM  REINITIALIZE  DATAFILE 
'3180  PRINT  #4H,GET(X,#2')  :REM  KEYWORD 

'3185  PRINT  #02,GET(X,#1) ,"  " ,GET(X , #1 ) , "  ";:REM  MONTH  DAY 

'3190  PRINT  #02,GET(X,#1) ,":",GET(X,#1) ,'':",GET(X,#1) ,"  ";:REM  HR:MM:SS 

'3200  PRINT  #02H,  GET(X , #1 )  :REM  COUNT 

'3205  FOR  L=1  TO  6:REM  T( 0 : 3 ) ,PITCH ,ROLL 
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'3206  PRINT  #04H,GET(X,#2) 

'3207  NEXT  LrPRINT 

'3210  FOR  M=1  TO  6:REM  2  VORT  SENSORS 
'3220  FOR  L=1  TO  4: REM  4  WORDS/LINE 

'3230  PRINT  #04H,GET(X,#2)  ";:REM  EACH  SENSOR  ON  A  LINE 
'3240  NEXT  L:PRINT:REM  FOUR  AXES/LINE 
'3250  NEXT  M 

'3262  PRINT  #3 ,GET(X,#2) , ”  ";:REM  COMPASS 
'3264  PRINT  #08H,GET(X,#4 ) , "  ";:REM  PRESSURE 
'3275  PRINT: REM  SPACE 
3280  PRINT  {0.,73}::BEM  TAKES  77 

3290  N=N+1 

3295  SLEEP  1>1:P=SDI (24): REM  SAMPLE  COUNTERS  AND  READ  PRESSURE 
3300  GOTO  3000 

8000  X=&H7300;REM  MULTIPLEXOR  LIST 


8010 

ASM 

X,DW 

&H0000;DW 

&H0101;DW 

&H0202;DW 

&H0303:REM 

ACM 

1/T1-T4 

8015 

ASM 

X,DW 

&H088A;DW 

&H098B;DW 

&H0AFF;DW 

&H0BFF:REM 

ACM 

1/PITCH, 

8020 

ASM 

X,DW 

&H20FF;DW 

&H21FF;DW 

&H22FF;DW 

&H23FF:REM 

ACM 

2 

8025 

ASM 

X,DW 

&H28FF;DW 

&H29FF;DW 

&H2AFF;DW 

&H2BFF 

8030 

ASM 

X,DW 

&H4000;DW 

&H4100;DW 

&H4200;DW 

&H4300:REM 

ACM 

3 

8035 

ASM 

X,DW 

&H4800;PW 

&H4900;DW 

J&H4A0p?pW 

fiH4BOO 

8040 

ASM 

x,im 

&H60b0;DW 

&H6100;DW 

<i(H6200  ;bW 

eiH630d:REM 

ACM 

4 

8045 

ASM 

X,DW 

&H6800;DW 

&H6900;DW 

&H6A00;DW 

&H6B00 

8050 

ASM 

X,DW 

&H8000;DW 

&H8100;DW 

&H8200;DW 

&H8300:REM 

ACM 

5 

8060 

ASM 

X,DW 

&H8800;DW 

&H8900;DW 

&H8A00;DW 

&H8B00 

8070 

ASM 

X,DW 

&HA000;DW 

&HA100;DW 

&HA200;DW 

&HA300:REM 

ACM 

6 

8080 

ASM 

X,DW 

&HA800;DW 

&HA900;DW 

&HAA00;DW 

&HAB00 

8090 

ASM 

X,DW 

&HFF00;DW 

&HFF00:REM  END  OF  T-IST 

8100  X=&H3800:A=X:REM  A/D  ROUTINE 

'  X  IS  POINTER  TO  AUX  ARRAY,  AB  IS  POINTER  TO  FLAGGED  CONVERSION  STORAGE 
8105  ASM  X,STD  &H7000:REM  POINTER  TO  DATAFILE 

8107  ASM  X,LDD  #&H6F38;STD  &H7002:REM  POINTER  TO  FLAGGED  CONVERSIONS 
8110  ASM  X,OIM  &H04,&H17;PSHX;LDX  #&H0600:REM  POWER  ON,  WAIT  5ms 
8120  D=X:ASM  X,DEX;BNE  D;PULX:REM  WAIT  (4* . 81us/LOOP) 

8150  ASM  X,CLR  &H11;LDAA  #&H10;STAA  &H10:REM  DISABLE  UART,  SET  SCI  FOR 
8160  ASM  X,OIM  8,&H11:REM  INTERNALLY  CLOCKED  8  BIT  DATA  RECEIVE  ENABLE 
'DUMMY  PULSE  SENT  OUT  TO  INITIALIZE 

8170  ASM  X,OIM  &HAB,&H157AIM  &HBF,&H15:REM  DUMMY  AUX  LOCATION 

8171  ASM  X,OIM  &H80,&H03;OIM  &H80,&H03;OIM  &H80,SeH03;OIM  &H80,&H03 
'START  TIMING  P27=l 

8172  ASM  X,AIM  &H7F, &H03 :REM  REMOVE  PULSE  P27=0 

8174  D=X:ASM  X,LDAB  #&H02;BITB  &H03;BEQ  D:REM  CHECK  P21  FOR  DONE 

8179  C=X:ASM  X,AIM  &H14,&H15:REM  ENTRY  POINT  FOR  AUX, MASK  PORT  5 

8180  ASM  X,LDAA  &H00,X;INX;INX:REM  TOP  OF  MUX  LIST, MOVE  TO  NEXT 
8185  ASM  X,ORAA  &H15;STAA  &H15:REM  PUT  MUX  ON  PORT  5 

8200  ASM  X,OIM  &H80  , &H03  ;OIM  &H80 ,&H03 ;OIM  &H80  , &H03 ;OIM  £cH80,&H03 
'START  TIMING  P27=l 

8210  ASM  X,AIM  &H7F, &H03 ;REM  REMOVE  PULSE  P27=0 

8220  D=X:ASM  X,LDAB  #&H02;BITB  &H03;BEQ  D;REM  CHECK  P21  FOR  DONE 

8240  ASM  X,LDAB  &H15;ANDB  #&H14;REM  MASKED  PORT  5  IN  B 

8250  ASM  X,LDAA  &H00  ,X;INX;INX:REM  GET  NEXT  MUX  WORD,  ADVANCE 

8255  ASM  X,ABA;STAA  &H15;REM  ADD  MUX  WORD  TO  MASKED  PORT  5,  STORE  P5 

8280  ASM  X,OIM  &H80 , &H03 ;OIM  &H80 , &H03 ;OIM  &H80 , &H03 ;OIM  &H80,&H03 
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START  TIMING  P27=l 

8290  ASM  X,AIM  &H7F,&H03:REM  REMOVE  PULSE  P27=0 

8400  ASM  X,CPX  0;CPX  0;CPX  0;CPX  0 

8402  ASM  X,CPX  0;CPX  0;CPX  0;CPX  0 

8403  ASM  X,CPX  0;CPX  0;CPX  0;CPX  0 

8404  ASM  X,CPX  0;CPX  0;CPX  0;CPX  0 

8410  ASM  X,LDAA  &H11;LDAB  &H12:REM  TRIGGER  A  READ 
8420  ASM  X,CPX  0;CPX  0;CPX  0;CPX  0:REM  16ECLOCK 

8430  ASM  X,LDAA  &H11;LDAB  &H12;REM  HIGH  BYTE  IN  B  AND  READ  NEXT 
8440  ASM  X,CPX  0;CPX  0;CPX  0;CPX  0:REM  16ECLOCK 

8450  ASM  X,LDAA  &H12;REM  LOW  BYTE  IN  A  DON'T  READ  NEXT 
'REM  FLIP  THE  BITS  OF  THE  WORDS 

8451  ASM  X,ROLA;RORB;ROLA;RORB;ROLA;RORB;ROLA;RORB 

8452  ASM  X,ROLA;RORB;ROLA;RORB;ROLA;RORB;ROLA;RORB;ROLA 

8480  ASM  X,TIM  &H20,&H03;BEQ  G:REM  TEST  BOTH  RECEIVED 
8488  ASM  X,XGDX;RORB;BCS  M:REM  BRANCH  IF  AUXILLIARY  (ODD  X) 

8490  ASM  X,ROLB;XGDX;CLRA;CLRB;BRA  G:REM  IF  NOT  RECEIVED , CLEAR 
8495  M=X:ASM  X,ROLB;XGDX 

'  PUT  AWAY  TO  §(50)  ARRAY  THE  FLAGGED  VALUES 

8510  G=X:ASM  X,PSHX;LDX  &H7002;STD  0,X;INX;INX;STX  &H7002;PULX 

8520  ASM  X,LDAB  )ff&HFF;EORB  &H00,X;BEQ  L;JMP  D:REM  CHECK  FOR  END  OF  LIST 

8550  L=X:ASM  X,LDAB  #&H02:REM  TEST  FOR  HOLD  PULSE 
8555  D=X;ASM  X,BITB  &H03;BEQ  D:REM  CHECK  P21 


8572  ASM  X,CPX  0;CPX  0;CPX  0;CPX  0 

8573  ASM  X,CPX  0;CPX  0;CPX  0;CPX  0 

8574  ASM  X,CPX  0;CPX  0;CPX  0;CPX  0 

8575  ASM  X,CPX  0;CPX  0;CPX  0;CPX  0 

8580  ASM  X,LDAA  &H11;  LDAB  &H12:REM  TRIGGER  A  READ 
8590  ASM  X,CPX  0;CPX  0;CPX  0;CPX  0:REM  16ECLOCK 

8600  ASM  X,LDAA  &H11;LDAB  &H12:REM  HIGH  BYTE  IN  B  AND  READ  NEXT 
8610  ASM  X,CPX  0;CPX  0;CPX  0;CPX  0:REM  16ECLOCK 

8620  ASM  X,LDAA  &H12:REM  LOW  BYTE  IN  A  DON'T  READ  NEXT 

'REM  FLIP  THE  BITS  OF  THE  WORDS  PUTS  HIGH  BYTE  IN  A,  LOW  BYTE  IN  B 

8621  ASM  X,ROLA;RORB;ROLA;RORB;ROLA;RORB;ROLA;RORB 

8622  ASM  X,ROLA;RORB;ROLA;RORB;ROLA;RORB;ROLA;RORB;ROLA 


8630  ASM  X,TIM  &H20 , &H03 ;BEQ  H:REM  TEST  BOTH  RECEIVED 
8640  ASM  X,XGDX;RORB;BCS  N:REM  BRANCH  IF  AUXILLIARY  (ODD  X) 

8645  ASM  X,ROLB;XGDX;CLRA;CLRB;BRA  H:REM  IF  NOT  RECEIVED , CLEAR 
8650  N=X:ASM  X,ROLB;XGDX 

'  FINISH  PUT  AWAY  TO  §(50)  ARRAY  THE  FLAGGED  VALUES 

8660  H=X:ASM  X,PSHX;LDX  &H7002;STD  0,X;INX;INX;STX  &H7002;PULX 

8670  ASM  X,XGDX;RORB;BCS  0:REM  BRANCH  IF  AUXILLIARY  (ODD  X) 

8680  ASM  X,ROLB;XGDX;OIM  &H40 , &H15 ;OIM  &H01 , &H17 :REM  PREPARE  FOR  AUX 
8690  ASM  X,LDX  #&H7301;LDD  &H7000;STD  &H7002;JMP  C:REM  DO  AUX 

'  EXIT  AFTER  DOING  ACM  AND  AUX 

8700  0=X:ASM  X,ROLB;XGDX 

8705  ASM  X,CLR  &H11:REM  DISABLE  SCI 

8710  ASM  X,LDAA  #5;STAA  &H10;REM  BAUD  RATE  FROM  TIMERl  (#5  GIVES  9600) 
8720  ASM  X,OIM  &H42 , &H1B:REM  ENABLE  TIMER2  INTS  AND  E/128 
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'8730  ASM  X,LDAA  #95;STAA  &H1C:REM  TIMER  COUNT  TO  GET  100  Hz 
8740  ASM  X,OIM  &H02 , &H11 :REM  02  RESTART  UART  NO  RECEIVE 
>  reconf  uart,  ENABLE  IT,  power  off  (PHYSICALLY  ENABLE  IT) 
8750  ASM  X,AIM  &H14 , &H15 :REM  PUT  MULTIPLEXORS  ON  PARK 
8760  ASM  X,AIM  &HFA,&H17:REM  P60 ,P62=0 , AUX  &  POWER  OFF 
8775  ASM  X,LDD  &H7002 
8780  ASM  X,RTS 


9100  B=X;REM  SUBTRACT  AND  TRANSFER  SUBROUTINE 

'  POINTERS:  m(7000)  =  DATA  ARRAY  PERMANENT  STORAGE.  PASSED  IN  AB 
'  in(7002)  =  FLAGGED  AXIS  PAIRS  STORAGE 

/  6F38  IS  @(50)  IN  RAM  MUST  BE  SAME  AS  A/D  ROUTINE 

'  in(  X  )  =  AUX  LIST  FOR  END  OF  LIST  CHECKING.  PASSED  IN  X 

9110  ASM  X,  STD  &H7000 

9115  ASM  X,  LDD  #&H6F38;  STD  &H7002 

9120  K=X;  ASM  X,PSHX;  LDX  &H7002;  LDD  0,X;  INX;INX:REM  GET  WORD 
9140  ASM  X,BEQ  F:REM  TEST  FLAG  ON  NORMAL  MEAS 

9150  ASM  X,SUBD  0,X:REM  DOUBLE  SUBTRACT 

9152  ASM  X,RORA;RORB;REM  SHIFT  RIGHT  WITH  CARRY  =  DIVIDE  BY  TWO 

9160  ASM  X,TST  0,X:REM  TEST  FLAG  ON  REVERSED  MEAS 
9170  ASM  X,BEQ  F 

9200  D=X  :  ASM  X,  INX;INX;STX  &H7002 ;PULX:REM  MOVE  ON  TO  NEXT  PAIR 

9203  REM  '  SAVE  RESULT  TO  DATAARRAY  AT  m(7000) 

9205  ASM  X,PSHX;LDX  &H7000;STD  0,X;INX;INX;STX  &H7000;PULX 

'  CHECK  FOR  END  OF  LIST.  AUX  LIST=  FOUR  BYTES  /  AXIS 
9240  ASM  X,INX;INX;INX;INX 

9250  ASM  X,LDAB  #&HFF:REM  CHECK  FOR  END  OF  LIST 
9255  ASM  X,EORB  &H00,X 

9260  ASM  X,BNE  K:REM  LOOP  IF  NOT  DONE 

9265  ASM  X,LDD  &H7000:REM  LOAD  DATAFILE  POINTER  FOR  RETURN 

9270  ASM  X, RTS: REM  EXIT 

9275  F=X 
'REM  FLAG  MISSED 
9280  ASM  X,LDD  #&H8000 
9290  ASM  X,BRA  D 

'  FIND  COMPASS  N+1  TRANSMIT  GAP 

9500  T=X:ASM  X,LDAA  &H03:REM  PORT  2 

9510  D=X:ASM  X,CMPA  &H03;BEQ  D:REM  FIND  EDGE 

9520  ASM  X,LDAA  &H03;LDX  #&H007B:REM  1ms  GAP 

9530  P=X:ASM  X,CMPA  &H03;BNE  D;DEX;BNE  P:REM  FIND  1ms  GAP 

9540  ASM  X,RTS 

9900  RETURN 


9999  END 
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